Non-Markovian master equation in strong-coupling regime 
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The time-convolutionless (TCL) non-Markovian master equation was generally thought to break down at 
finite time due to its singularity and fail to produce the asymptotic behavior in strong coupling regime. However, 
in this paper, we show that the singularity is not an obstacle for validity of the TCL master equation. Further, we 
propose a multiscale perturbative method valid for solving the TCL master equation in strong coupling regime, 
though the ordinary perturbative method invalidates therein. 
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I. INTRODUCTION 

The study of non-Markovian quantum open system attracts 
increasing attention nowadays. There are two reasons for 
the interest. On the one hand, the popular Markovian 
approximation which neglects the memory effects of the 
environment is not sufficient for the recent progress in many 
fields, such as quantum information processing JW, quantum 
optics condensed matterphysics [g], chemical 

physics [7], and even life science [8]. On the other hand, there 
are still many open questions for the theory of non-Markovian 
quantum open system. 

Several methods are proposed to study the non-Markovian 
open quantum systems J2, IH-tHl], among which the non- 
Markovian master equation is quite promising. Projection 
operator techniques provide systematic framework to derive 
master equations. Different projection operator techniques 
give different kinds of non-Markovian master equations. Two 
kinds of non-Markovian master equations, the Nakajima- 
Zwanzig master equation |9] and the time-convolutionless 
(TCL) master equation 111 Oil , are widely used. Since the 
Nakajima-Zwanzig master equation is an integro-differential 
equation, the TCL master equaiton which is a time-local 
first order differential equaiton is much easier for numerical 
solution. Besides the methods by extending the Hilbert space 
EH H, a new numerical method called non-Markovian 
quantum jump |15] and its modified scheme [16] were 
proposed recently. 

However, in strong-coupling regime, it was generally 
thought that there are two severe problems for the application 
of the TCL master equation. One is that the TCL master 
equation breaks down at finite time in strong-coupling regime 
due to the singularity, thus fails to produce the asymptotic 
behavior Jj, fl711 . Another is that the ordinary perturbative 
method fails to produce the correct behavior 0]. 



In this Letter, we study the dynamics of a two-level system 
interacting with a structured environment in strong-coupling 
regime. Our result shows that the singularity at finite time 
is not an obstacle for the TCL master equaiton to produce 
the correct asymptotic behavior. Moreover, we introduce a 
multiscale perturbative method which produces the correct 
behavior, though the ordinary perturbative method fails, in 
strong-coupling regime. 



II. THEORETICAL FRAMEWORK 

Under rotating wave approximation, the total Hamiltonian 
of the two-level system with the bosonic reservoir in zero- 
temperature is given by H = Hs + He + Hj = Hq + 
Ht , with H s = w cr + o-_ , H E = Xk ^k a \ a k > H i = 
Tik (gkC+a/c + g* k a k <T-) , (H = 1). Here <x ± and ojq are the 
inversion operators and transition frequency of the two-level 
system, respectively, a' k , the creation and annihilation 
operators of the field modes of the reservoir with frequency 
a>k, and gk, g* k the coupling strength between the two-level 
system and the kth field mode of the reservoir. Since [H, N] = 
0, where N = cr+cr- + 2* ctiak, for an initial state of the form 
| (A(0)> = (c g o\ g) + c e o\ e))\ 0>£, the time evolution of the 
total system is confined to the subspace spanned by the bases 
{IS>|0> £ ,|e>|0> £ ,ls>|l t > £ }as 

I m) = Cgo\ g)\ 0)e + c e (t)\ e)\ Q) E + J] Q (f)| g)\ We , (D 



where | \k)E is the state of the reservoir with only one exciton 
in the fcth mode. 

According to the Schrodinger equation in the interaction 
picture with Hq = H$ + He, one can obtain an integro- 
differential equation for the amplitude c e (f) as 
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FIG. 1 : (Color online) Time evolution of the population in the upper 
level for different initial states (a) | (f>\) = \ e) and (b) | $2) = (I e) + 
I g»/ V2 , (t = At), both with A = 1, y = 10, for the cases of (i) the 
numerical solution of TCL master equation for initial state | </>] > = | e) 
(asterisk), (ii) the exact analytical solution for initial state | 0i) = | e) 
(solid curve), (iii) the numerical solution of TCL master equation for 
initial state | <p 2 ) = (| e) + | g))/ V2 (diamond), and (iv) the exact 
analytical solution for initial state | <p 2 ) = (I e) + | g))/ V2 (dashed 
curve). 



where the correlation function f(t - t) takes the form 

fit - r) = V I gk \V<*>-<*^ = J{a> k y {n)i, - T) da> k , 

where J{a)k) is the spectral density function of the reservoir. 
Thus, the TCL master equation takes the following form (2, 

\vh\m 



dt 



Ps(t) 



--S(.t)[o-+cr-,Ps(t)] 

+y(t)[a-p s (t)o-+ - i{cr + o-_,p,(*)}] , (3) 



with Lamb shift and decay rate given by 

c e {t) 



S(t) = -21m 



c e (t) 



, y(t) = -2Re 



C e (t) 



Ce(t) 



(4) 



two-level system interacts with the central frequency of the 
reservoir resonantly, loq = co c J2,[T3l. In the following, we 
restudy this problem for this model. If a> c » A, such as in 
an optical cavity, co c can be extended to infinity. Then the 
correlation function f(t - r) is given by 



ft* ^ - 
fit - t) = — e 



At-A 



(5) 



Substituting Eq. (0 into (O and using Laplace approach, one 
obtains J2[] 



c e (t) = c e0 e 



-At/2 
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+ ^sin 

















(6) 



where T = y/2yoA - A 2 . For the strong-coupling case (yo > 
A/2), F > 0, and c e (t) is an oscillating function with discrete 
zeros at t n = 2[(n+ l)n - arctan(T/A)]/r, (n = 0,1,2,--)- 
Substituting Eq. © into (0]), one obtains the exact expressions 
for S(t) and y(t) im 



5(0 = 0, y(t) 



2y tan(rf/2) 
1 +i/rtan(n/2)' 



(7) 



A. Validity of TCL master equation 



Therefore, we see that y(t) diverges at these points t„. 

Previously, it was generally thought that the singularity is 
an obstacle for validity of the TCL master equation 1121 UJn . 
On one hand, for TCL master equation, since "the evolution 
of the reduced density matrix only depends on the actual 
value of p s (t) and on the TCL generator" [2] and the density 
matrices coincide at t — to for different initial states, the 
evolution of the density matrices after to should be the same 
for different initial states. On the other hand, the exact 
analytical solution ||2|] for the problem shows that the density 
matrices for t > to differs for different initial states. That 
means the solution of TCL master equation does not agree 
with the exact analytical solution for t > to- So, it was thought 
that, "...a time-convolutionless form of the equation of motion 
which is local in time ceases to exist for t > to..!' 10] or 
"...the time-convolutionless generator breaks down at finite 
time in the strong coupling regime, thus failing to reproduce 
the asymptotic behavior..." 11711 . 

In our opinion, the singularity is not an obstacle for validity 
of the TCL master equation. By using the method in llrjll . 
we solve the TCL master equation Eq. ([3]) numerically. In 
the simulation, the decay rate (Eq. (0) with singularity is 
used. From Fig. Q] we find that the numerical solution of TCL 
master equation agrees with the exact analytical solution very 
well for t > to- That means, even though the TCL master 
equation has a singular point at t = to in strong-coupling 
regime, it still reproduces the correct dynamics when t > to- 
Actually, since t = to is a singulary point of the TCL master 
equation, the dynamics around t — to cannot be explained by 
the theory of the first-order ordinary differential equation at an 
ordinary point lfl9ll . 



One problem for the TCL master equation is that it was 
thought to break down at finite time in strong-coupling regime 
due to the singularity. For example, the problem occurs in 
the case that the spectral density of the reservoir takes the 
Lorentzian form J{<x>k) = yoA 2 /2jt[(cl)i 1 - a> c ) 2 + A 2 ] and the 



B. Multiscale perturbative expansion 

Another crucial problem for TCL master equation is that 
the ordinary perturbative expansion fails in strong-coupling 
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FIG. 2: (Color online) Time evolution of the population in the upper 
level for initial state | e) by multiscale perturbative method and 
ordinary perturbative method (t = At), with A = 1, 70 = 10, for 
the cases of (i) the exact solution (solid curve), (ii) the numerical 
solution of exact TCL master equation (cross), (iii) the first order 
multiscale method (diamond), (iv) the second order multicale method 
(asterisk), (v) the second order ordinary perturbative method (dash- 
dotted curve), (vi) the fourth order ordinary perturbative method 
(dotted curve). 



regime [2]. The reason is that the ordinary perturbative 
expansion corresponds basically to a Taylor expansion of y(t) 
in powers of jq. In fact, this method treats the dynamics only 
in one time scale. For the model considered above in strong- 
coupling regime, there are two time scales, which correspond 
to the decaying and the oscillating behaviors, respectively. 
The ordinary perturbative expansion only considers the time 
scale of the decaying behavior, so the oscillating behavior 
disappears in the perturbative solution (see Fig. 2). 

Since the failure of the ordinary perturbative expansion 
originates from ignoring multiscales of the dynamics, we 
introduce a multiscale perturbative expansion 12011 to treat the 
strong-coupling case. 

According to Eq. by giving a multiscale perturbative 
expansion of c e {t), one can get the multiscale perturbative 
expansions of y(t) and S (t). From Eqs. (|2]) and (0, one obtains 



c e (t) + Ac e (t) + ?j-c e (t) = . 



(8) 



By introducing dimensionless parameters T = yot and 
e-A/yo, where s <K 1 for strong-coupling regime, Eq. <j8j 
reads 



with initial conditions c e (0) — c e o and c e (0) = 0. There 
are two kinds of behaviors of the dynamics, corresponding 
to two different time scales. The time scale of the decaying 
behavior relates to At = sT, and that of the oscillating 
behavior relates to a complicated function of s. Therefore, 
two different time scales, t\ = (e'^aj + s 3 ^ 2 a 2 + and 
t 2 = sT, are introduced, where a,'s are unknown parameters 
to be determined. Expanding c e {t\, t 2 ) in powers of s 



c e (tut 2 ) = cf\t u t 2 ) + s l,z c^(tut 2 ) + ecf>(tut 2 ) + ... , (10) 



,1/2,(1) 



and substituting Eq. ( TTOb into (0 and the initial conditions, 



one obtains the equations for as follows 



To the order of s, one gets the equations for cf^ as 



fl 2iL c (0) + i c (0) = , 

l dt\ 2 e 



cf ) (0,0) = c e0 , £-4°>(0,0) = 0. 

Ot\ 



To the order of s 3 ^ 2 , one gets the equations for c e as 



d 2 (0) d „ 

Ot 1 Ot 2 Ot\ 



dt\ 2 
0^(0,0) = 0, 



To the order of e 2 , one gets the equations for Cg as 



ft 1 1 ff- 8 2 

2_£_ (2 A (2 _ r. (1) r. (0 

ai dtf e + 2 Ce ~ 2ai d h dt 2 C * 2aiai dt\ Ce 

_^l c «»)_ fl ,A c o)_ JL C «>) 
dt\ e l dt x Ce dt 2 Ce ' 



c< 2) (0,0) = 0, 



a\ t-4 2) (0, 0) + a 2 ^-cf\0, 0) + — c^fO, 0) = . 



d_ 



d_ 

'dt[ 



By a routine multiscale analysis of the above equations 
ll20ll . one obtains the solutions of By substituting the 
perturbative solution of c e into Eq. (0), we can get the 
corresponding Lamb shifts and decay rates. The first order 
solution is 

S(t) = 0, y(t) = e + V2^tan( ^T/2), 



d d e 

dT-° e + E df Ce + 2° e 



(9) 



and the second order solution is 



4 



^f^{e 3/2 cos[ V2e(l - e/4)T/2] + V2(4 + e) sin[ V2e(l - s/4)7//2]} 
4cos[ V7s(l - e/4)T/2] + 2 V2^sin[ V2e(l - e/4)T/2] 



By solving the corresponding master equations using the 
method in itlol . we study the dynamics of the population 
in the upper level. From Fig. [2] we can see that, unlike 
the solutions of ordinary perturbative method where the 
oscillating behavior is missing, the decaying and oscillating 
behaviors are both included by multiscale perturbative 
solutions. The perturbative solution up to the second order 
fits very well with the exact solution. 

For the exact solution, Eq. (0, the first singular time of 
y(f) is fo = 2 arccos(- Vs/2)/[yo V(2 — s)s]. For the first 
and second order approximations of fo by the multiscale 
perturbative expansion, the first singular time of y(t) are 
7T/[y ^/2e] and 4 V2arccos[- Ve/(2 ; + s)]/[yo(4 - e)Vi], 
respectively. Detailed analysis shows that relative errors of 
the first and second order approximations are in the orders of 
s 1 / 2 and e 3 ^ 2 , respectively. 

III. CONCLUSION 

To summarize, we study the time-convolutionless non- 
Markovian master equation in strong-coupling regime. For 
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